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Resistance and reactance of diathermy patients to 
alternating currents of varying frequencies up to and 
including radiofrequencies have been measured with a 
high frequency Wheatstone bridge. At frequencies greater 
than one million cycles per second the resistance becomes 
constant and the reactance and phase angle decrease to 
zero. Under these conditions the human body behaves as 


a pure resistance and the heat production in calories can 
be computed from the formula 0.24 VJt where V is the root- 
mean-square voltage, J the root-mean-square current and 
t the duration time in seconds of the diathermy treatment. 
High frequency resistances of various parts of the body 
which are usually treated with diathermy are given. 


N medical diathermy a high frequency alter- 

nating current passes between electrodes 
applied to the skin of the patient and heat is 
generated in the tissues traversed by the current. 
The frequency of the alternating current is 
sufficiently great to prevent electrical stimulation 
of the neuromuscular system. In physiological 
problems related to temperature-control and 
water-balance it is desired to measure the total 
amount of heat imparted to an animal or to a 
patient during a diathermy treatment. It is 
obvious that calorimetric measurements of this 
quantity would be extremely involved. A simple 
method would be to measure the high frequency 
current and high frequency voltage and, if the 
phase angle were known, it would be possible to 
compute the heat production from these quan- 
tities. The present paper is concerned with a 
measurement of the electrical resistance and 
reactance of the patient during a diathermy 
treatment at frequencies up to and including 
diathermy frequencies. Knowing the resistance 
and reactance at various frequencies it is possible 
to compute the phase angle. 

There is evidence that the phase angle of 


biological tissues decreases to a low value at 
radiofrequencies. Philippson! found by measuring 
impedance of tissues at various frequencies that 
the impedance became constant for values of the 
frequency greater. than one million cycles per 
second. Fricke? measured the resistance and 
capacity of red cell suspensions at various fre- 
quencies and found that the resistance decreased 
to a low constant value between one and ten 
million cycles per second. This constancy of the 
resistance at high frequency indicates that the 
phase angle is small at high frequency. McClen- 
don*: * observed a large decrease of resistance of 
red cell suspension as the frequency of the a.c. 
current used to measure resistance increased. 
Hemingway’ in calorimetric determinations of 
the heat production by diathermy currents in 
vegetable tissue found that the phase angle of 
the diathermy current was very small. 
Measurements of resistance and capacity were 


! Philippson, Comptes Rendus Soc. Biol. 83, 287 (1920). 
? Fricke, J. Gen. Physiol. 9, 153 (1925). 

3 McClendon, J. Biol. Chem. 69, 733 (1926). 

4 McClendon, Am. J. Physiol. 82, 525 (1927). 

5 Hemingway, Radiology 17, 136 (1931). 
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made on the high frequency Wheatstone bridge 
described in a previous paper by Hemingway and 
McClendon.* Additional coils and capacities have 
been inserted into the oscillating circuits and the 
detector circuit, in order to measure resistances 
and reactances with alternating currents at fre- 
quencies from the audiorange of frequencies 
upwards. In order to eliminate the use of an 
extremely high capacity at lower frequencies, 
the arrangement of the balancing arm of the 
bridge was a variable high frequency resistor in 
parallel with a variable air-mica condenser. The 
corresponding series resistances and capacities 
were computed from the observed parallel 
resistances and capacities by the well-known 
formulas 


Series resistance = R/(1+S*) ohms, 
Series reactance = RS/(1+.S*) ohms, 


where R is the observed parallel resistance, C is 
the observed parallel capacity, f is the frequency 
of the alternating current and S=27fRC. 
Measurements of resistance and capacity were 
made with alternating currents of frequencies 
from 1000 cycles per second to 1.5X10* cycles 
per second. The resulting series resistances and 
reactances of the tissues with electrodes at 
various parts of the body are given in Fig. 1. It 
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Fic. 1. Variation of diathermy resistance and reactance 
with frequency. 


* Hemingway and McClendon, Physics 2, 396 (1932). 


is to be observed that, as the frequency of the 
alternating current increased, the reactance 
decreased and approached zero in the range of 
frequencies which include diathermy alternating 
currents. In this range the resistance approached 
a constant value. Hence, in this diathermy range 
of frequencies the phase angle approaches zero. 

Only three representative curves showing the 
variation of reactance with resistance are given. 
Others have a similar shape but have been 
omitted to prevent confusion because of too 
many curves on a single graph. The knee-ankle 
resistance and reactance values were obtained 
from a cuff electrode above the knee and an 
electrode below the sole of the foot, which was 
immersed in salt solution. For the knee, cuff 
electrodes above and below the knee were used. 
For the shoulder, two 9X 16 cm? electrodes were 
used, one anterior, the other posterior. The fol- 
lowing table contains values of the high fre- 
quency resistance of patients measured with 
alternating current of 1.5 million cycles per 
second, together with the type and size of the 
diathermy electrode applied. 


TABLE I, 

Region Resistance 

treated Electrode (ohms) 
Shoulder Metal 9X16 cm? 37 
Shoulder Metal 9X16 cm? 28 
Shoulder Metal 9X16 cm? back) 40 

or Metal 9X16 cm? (back) \ 
Sacro-iliac {Metal 17X17 cm? (front)} 34 
Knee Pad 9X 16 cm? -40 
Knee Pad 9X 16 cm? 50 
Ankle (calf to sole) Cuff and pad 150 


DISCUSSION OF RESULTS 


The curves of reactance and resistance vs. 
frequency show that as the diathermy range of 
frequencies is approached the phase angle, i.e., 
arctan (X/R) where X is the series reactance 
and R is the resistance, becomes zero. The heat 
production in calories is given by the expression 
0.24 VJItcos 6 where V is the root-mean-square 
voltage and J the root-mean-square current and 
t the time in seconds. Since @ is zero, or very 
close to zero, the value of the heat production in 
calories during a time of t seconds becomes 0.24 
VIt. This shows that the human body behaves 
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to diathermy current as a pure resistance and the 
heat production can be computed from the 
observed voltage and current values. 

Diathermy machines are equipped with a 
meter which indicates the root-mean-square 
current and in machines examined and tested it 
has been found that the meters are quite accu- 
rate. The measurement of the root-mean-square 
high frequency voltage presents certain dif- 
ficulties in regard to balanced shielding and 
elimination of stray currents. The Weston 
Instrument Company have recently designed 
two voltmeters on the thermocouple principle 


for the purpose of measuring high frequency 
diathermy voltage, which have proved satis- 
factory. Full scale deflections of these instru- 
ments are 50 and 125 volts and measurement of 
heat production in animal tissues during dia- 
thermy treatments is being carried on with the 
use of these instruments. 

One of these instruments was purchased with 
a grant from the American Medical Association 
and we wish to take this occasion to thank the 
Association and Mr. H. A. Carter, Secretary of 
the Council of Physical Therapy, for this 
assistance. 
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Note on Dipole Radiation Theory 


W. Howarp Wise, American Telephone and Telegraph Company, New York 
(Received June 17, 1933) 


By fitting certain general asymptotic series solutions the upper infinite hemisphere, recurrence relations between 
of the wave equation to known expressions for wave the successive terms of the asymptotic series for some of 
function components at the surface of the ground and on _ the dipole wave function components are obtained. 


INTRODUCTION 


AST attempts to calculate the radiation field of a dipole on or near a flat homogeneous earth 

have been chiefly confined to efforts at direct evaluation of the Sommerfeld and von Hoer- 
schelmann infinite integral representations of dipole wave functions. Because of the great com- 
plexity of these integrals this direct attack has failed to yield an explicit recurrence relation for any 
of the wave function component series in the air. However, in the case of the vertical dipole of zero 
height, centered on the interface between air and ground, rigorous series expansions, both convergent 
and asymptotic, which hold at the surface of the ground, have been derived for each of the wave 
function components! and the first term of the asymptotic expansion for the sky wave from any 
dipole is known.? In these cases where we have rigorous expressions for wave function components 
on the boundaries of the air (the surface of the flat earth and the upper infinite hemisphere) the 
problem may be regarded as one of finding the solutions of the wave equation in air which reduce to 
these expressions on the boundaries. In this paper recurrence relations for two of the wave function 
component series are obtained by fitting appropriate asymptotic series solutions of the wave equation 
to two of the known expressions for these wave function components on the boundaries of the air. 
The solutions of the wave equation assumed for this purpose appear to be new. 


A RECURRENCE RELATION FOR Q; 


In the air, the wave function for a dipole in air is a solution of the wave equation V?W+k?W=0. 
The primary wave function is the simple solution exp (tkir)/r. The whole wave function is 
u+exp (ikir)/r, where u is the wave function for the currents in the ground. Then V*u+k,u=0. 

Sommerfeld has shown that u may be conveniently separated into three component parts, 
u=P+Q:1+(Qs2, of which P is the Zenneck wave or ground wave which is completely known and so 
not considered here, Q; is a spherical wave and Q, is a highly attenuated, and therefore for most 
purposes relatively insignificant, wave of hybrid character. 

In polar coordinates (see Fig. 1) V?u+ku=0 is 

1 Ou 1 ou 1 
— — rf — sin 0 —_+ =0. 
Ore Ore re sin 8 00 00 re? sin? 6 


Writing ikir2= and cos @=c this takes the somewhat simpler form 


1 ou 10 Ou 1 
— — & —+— — (1-c?) —+—_- —_-_ 0. (1) 
0a dc dc @(1-—c) 
1W. H. Wise, Proc. 1. R. E. 19, 1684 (1931). 2? W. H. Wise, Bell Sys. Tech. J. 8, 662 (1929). 
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Pp 
Z 
+ 
Fic. 1. 
Assuming a solution of the form - 
= 1k, cos m(¢—¢o) (2) 
one gets 
o m?* 
tk, cos m(¢—¢o) (1 —c*) gmn'’(c) = 2cgZmn' (Cc) n(n 1) — 2ngmin41)(C) =0, 
n=1 @t 2 
whence 
@-(— ) (3) 
m(n mn\C) ——Zmn \C mn \C). 


This is as would be suspected; since the first term of the asymptotic Q:-component series (2) is Q: 
on the sphere at infinity, a knowledge of it alone serves to completely determine Q,. 

Parenthetically, it may be remarked that the series (2) must be Q: and not the complete solution 
because P and Q: are exponentially attenuated and a series of this character cannot represent an 
exponentially attenuated function. 

In the case of the vertical dipole we know by previous work,’ and also by reciprocity, that the 
first term of Q; is R; exp ikire/re, where R; is the coefficient of reflection for a plane, electromagnetic 
wave polarized with the electric vector in the plane of incidence. Hence, in this case, m is zero and 


c+nr(1—7(1 


gu(c)=Ri= 
where 7 T=hki/ke and 
= enw? +4rwopiw =k; in the air (ki ~27/d) 
=k,? in the ground. 
In the case of the horizontal dipole we know by previous work? that the first term of Q; is the vector 
—2 cos etkire 


aX +bxX0+c¢X 
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where the a, 6 and c preceding the multiplication signs are unit vectors pointing in the x, y and z 
directions. The earth permeability has been assumed to be unity. The dipole points in the positive 
x-direction. Clearly, the x-component of Q; is defined by m =0 and 


m+(1—7°(1—<*))! 


g0(c) = 


while the z-component of Q, is defined by m=1, ¢9=0 and 


gu(c) = 


A TRANSFORMATION OF THE SERIES (2) 


Since the solution of the wave equation in polar coordinates as usually derived is a series of 
products of Hankel functions and associated Legendre functions, it is natural to seek to identify 
the series (2) with the ordinary formal solution. The ordinary type of solution appropriate for the 
representation of Q, is 


=ik; cos m(¢—¢o) ) H® (4) 


2kire 


1 
= tk, cos m(¢— go) Ao—Po™(c) +Ai1— (1-—) Pino: 
L @ @ @ 


n(n+1) 


€ 
= ik, cos m(¢—¢p) ill A,P.™(c)+— A,P,™(c) 
 (n—1)n(n+1)(n+2) 
A,P,»™(c)+:° ‘| 
2 2?X2! 


Now 


s—1 m? cd @& 
2 (1—c*)2s s dc 2s dc 
1 d? d m? 
=— —c*) —P,"™(c) —2c — {s(6- 1)- | 
dc? dc 


2s 


1 (n+1—s)(n+s) 
=— [s(s—1) —n(n+1) JP,™(c) = — Pr™(c). (5) 
2s 2s 


It follows that, if gmi(c)=>LA,»P,”(c) then, 
0 


n(n+1) wo (n—1)n(n+1)(n+2) 
&m2(c) = — AnPn™(c); &ma(c)= 
1 2x1! 2! 
Hence the two solutions are identical if the constants A, are defined by 
2n+1 (n—m)! 
(6) 
2 (n+m)!/_, 
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gZoi(c) is only defined for positive values of c; for negative values of c it is at our disposal. In the case 
of the vertical dipole the integration can be effected without much trouble. Since P,(c) is a poly- 
nomial in ¢ the chief problem is the evaluation of the integral 


1 
(7) 
0 
The integration is straightforward but a bit long. 
If goi(—c) = —2—goi(c) and 7=1/4 then 
A,= 42.1137, A;= —0.7795, A; =+0.4274, 


A; = —0.2689, Ay=+0.1829, Ay = —0.1311, 
So, at least for this particular choice of go:(—c), the series (4) is rather slowly convergent. 
An unfortunate property of the goi(—c) chosen is that while it makes go;(c) and its first derivative 


continuous at c=0 it also makes the second derivative and so also go2(c), discontinuous at c=0. 
This is a difficulty which it is not easy to avoid. 


A RECURRENCE RELATION FOR Q2 


It has already been found that at the surface of the ground the Q2 component of the wave function 
for a vertical dipole of zero height is! 


—~+-—+4+-— 


2 


tkir sD, De Ds; 
(8) 


3 


where n=ikep, Di=(1+7?), and 
We seek a solution of the wave equation which will reduce to (10) at the surface of the ground. 
To this end we assume the asymptotic solution form (9), 


tkit D, 


where (= k2z, substitute (9) into the wave equation in cylindrical coordinates, equate to zero the 
coefficient of each power of 1/n, solve the resulting sequence of differential equation in the ¢,(¢) 
and apply the boundary conditions ¢,(0)=1, ¢,(*) =0. 

The differential equations to be solved are 


— (1-7?) ($) = 
— (1-77) J = —3DidilS), 


It is easy to solve a few of these equations but, after the third, the labor required soon reaches 
prohibitive proportions. The solutions of the first two equations are 


=(14—. 
Dz 2(1—7?)! 


358 W. HOWARD WISE 


If a Dodo(f) term is added to (9) it turns out that ¢o(¢) =0. 
The Q2 just derived is the Q2 for the integral 


1+7? 1 
f Hy™ 
2 


where /= (\*—k,*)! and m= (\*—k,*)!, the wave function for a vertical dipole of zero height. 
If the vertical dipole is of finite height its wave function is 


etkir eikir2 1 
-—+ f Ho "dd. 
r -ol+r?m 


So, if the vertical dipole is of finite height its Q2 is (9) with a factor 2/(1+ 7°) added and z replaced 
by z+a. 
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The Effect of Inductance on the Intermittent Glow Discharge 


Winston E. Kock,* University of Cincinnati 
(Received May 20, 1933) 


The insertion of an inductance in the condenser arm of 
an intermittent glow discharge circuit was found to affect 
the characteristics of the oscillations. The voltage-fre- 
quency curve was noticeably altered, the wave form of 
the voltage across the condenser became sinusoidal and 


other new wave forms were introduced. Within a certain 
range, the frequency of oscillation was practically inde- 
pendent of the voltage, resistance or tube characteristics. 
Various applications are suggested. 


INTRODUCTORY 


HE intermittent glow discharge, first 

noticed by Gassoit! and later discussed by 
Hittdorf,? is obtained by placing a condenser in 
parallel with a discharge tube and a resistance in 
series with it and a voltage supply (Fig. 1). The 
approximate wave form of the voltage across the 
condenser is shown in Fig. 2.* The frequency of 
intermittence is approximately given by* 


f= (RC log. (V—Vo)/(V—V1))>, 


where V is the applied potential, Vo is the igni- 
tion potential and V, the extinction potential of 
the discharge tube. 

The inclusion of an inductance in the con- 
denser arm has been tried by several workers in 
this field but the effects with which this report 
is concerned were not recorded. 

Herweg? first inserted the inductance to check 
the predictions of Kaufmann® regarding the 
stability of the discharge. He utilized the circuit 
shown in Fig. 3 and, starting with a low value of 
R, he gradually increased this resistance until 
oscillations were produced. Over a very small 
range of R he obtained oscillations of the 


* Special Fellow in the Department of Physics, Uni- 
versity of Cincinnati. 

1 Gassoit, Proc. Roy. Soc. (London) Al2, 329 (1862). 

? Hittdorf, Ann. d. Physik 7, 566 ff. (1879). 

*’Oschwald and Tarrant, Proc. London Phys. Soc. 36, 
262 (1924). 

4See Pearson and Anson, Proc. London Phys. Soc. 34, 
204 (1922); also R. Anthouard, Comptes Rendus 185, 
188 (1927), 

5 Herweg, Phys. Zeits. 13, 633 (1912). 

® Kaufmann, Ann. d. Physik 307, 175 (1900). 


DISCHARGE TUBE 


R 


D.C. 
VOLTAGE SUPPLY 


Fic. 1. Circuit diagram for the intermittent glow discharge. 


first kind, i.e., oscillations wherein the current 
through the discharge tube never ceased com- 
pletely; and as R was further increased, he 
obtained intermittence oscillations of the second 
kind, oscillations wherein the current through 
the discharge tube was truly intermittent. 


VOLTAGE 


TIME 


Fic. 2. Wave form of the voltage across the condenser in 
the intermittent glow discharge. 
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DISCHARGE TUBE 


R 


VOLTAGE SUPPLY 


Fic. 3. Circuit diagram of the inductive glow 
discharge oscillator. 


Dallenbach’ later pointed out that Kaufmann’s 
formula was not applicable to Herweg’s circuit 
unless certain assumptions were made. In 1924, 
Oschwald and Tarrant’ noticed that at high 
frequencies of intermittence, the inclusion of an 
inductance in the condenser arm had the effect 
of maintaining oscillations above the limiting 
frequency otherwise obtainable. Clarkson* re- 
ported that inductance in the condenser arm 


a: 
800 A 
600P 
e 
z 
= x x 
400) 
% 
> 
200 
80 100 120 


SUPPLY VOLTAGE 


Fic. 4. Frequency-voltage curves; A, the intermittent 
glow discharge (no inductance, circuit of Fig. 1, C=0.001 
mf, R=1 megohm); B, the inductive glow discharge 
oscillator (circuit of Fig. 3, L=150 henries, C=0.001 mf, 
R=1 megohm). 


7 Dallenbach, Phys. Zeits. 27, 101 (1926). 
§ Clarkson, Phil. Mag. 4, 849 (1927). 


facilitated the production of a series of flashes 
of decreasing intensity. However, no mention 
has been made of the effect of inductance on the 
frequency-voltage curves or on the wave forms 
of the various voltages throughout the circuit; 
these effects will now be discussed. 


EXPERIMENTAL 


In this investigation the inductances were of 
the order of several hundred henries; for this 
purpose, the secondaries of small radio trans- 
formers (audiofrequency) were admirably suited. 
Commercial cold-cathode one-watt neon glow 
lamps were used for the glow discharge tube and 
batteries were employed for the voltage supply 


VOLTAGE ACROSS C 


VOLTAGE ACROSS R 


VOLTAGE ACROSS L 


Fic. 5. Wave forms of the inductive glow discharge 
oscillator when operated in the region X —X, Fig. 4. 


when rectified alternating current was not 
suitable. Measurements of frequency were made 
with a Western Electric oscillator and checked 
with tuning forks; wave forms were observed on 
a cathode-ray oscillograph and recorded on a 
Duddell oscillograph. 

Frequency-voltage curves are shown in Fig. 4. 
These curves apply to the circuit of Fig. 3; curve 
A portrays the characteristic with the inductance 
short-circuited, whereas curve B shows the 
characteristic with the inductance operative. For 
curve B the limiting frequency of intermittence 
is that frequency at which the oscillating circuit 


« 
A 
= 
id 


INTERMITTENT GLOW DISCHARGE 


resonates and corresponds to the frequency of 
Herweg’s® oscillations of the first kind, i.e., 
where C’ includes the inter- 
electrode capacity of the tube and L’ includes the 
effective inductance of the tube. The frequency 
is seen to be practically independent of the volt- 
age over a large region (portion X—X); if the 
frequency is plotted against values of the series 
resistance, the same flat characteristic is ob- 
tained. It is also evident that ordinary changes 
in the tube constants will not appreciably affect 
L’ and C’, hence, for the region X—X, the 
frequency is little affected by variations in 
voltage, resistance or tube constants. 

Various wave forms are shown in Fig. 5. 
Comparison with Fig. 2 indicates the change 
which the voltage across the condenser undergoes 
when the inductance is inserted. The new wave 
is very nearly sinusoidal if the tube is operated 
in the region X—X, Fig. 4.. However, the voltage 
across the series resistance is not greatly affected, 
because of the absence of inductance effects in 
the outer branch circuit. 
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SUGGESTED APPLICATIONS 


For a.c. bridge measurements and other 
experiments where harmonics in the voltage 
wave are objectionable, the inductive glow dis- 
charge oscillator offers a convenient and simple 
sine wave generator of practically constant 
frequency. Inasmuch as the voltage swing on the 
condenser is of the order of five to ten volts, only 
one stage of amplification would be required to 
provide a fair power output. Variations in fre- 
quency can be introduced by a variable con- 
denser, thus affording a sine wave oscillator of 
large frequency range. 

The inductive glow discharge oscillator is 
being utilized by the author for the electrical 
production of music. The remarkable pitch 
stability over the region X—X, Fig. 4, makes this 
circuit an excellent source of vibrations for an 
electro-musical instrument, wherein a separate 
discharge tube is employed for each note of the 
scale. A description of this instrument will appear 
in a later article. 
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The Temperature Distribution of a Gas Flowing Through a Furnace 


K. F. HERZFELD, The Johns Hopkins University, Baltimore, Maryland 
(Received July 27, 1933) 


In recent experiments the heat of chemical activation’ 


(formation of radicals) has been determined by measuring 
temperature coefficients of reaction rates of gases flowing 
through a furnace. To justify this procedure, an investi- 
gation is necessary to find out to what extent temperature 
equilibrium is reached. In this paper it is shown that the 
flow is laminar, and that two cases must be distinguished. 


Temperature equilibrium is reached throughout the flow 
when the pressure drop is less than 0.003 mm Hg per cm. 
For pressure drops about ten times as high, only a thin 
sheet close to the wall reaches temperature equilibrium, 
but this fact affects only the total number of activated 
molecules reaching the wall, not the temperature coeffi- 
cients. 


N recent experiments, the behavior of gases 

flowing through a cylindrical furnace at low 
pressure and high speed has been studied. The 
number of radicals produced in the furnace has 
been used as a basis on which to calculate the 
heat of activation. It is, of course, of primary 
importance to determine whether temperature 
equilibrium has really been attained. 


A. THe FLow 


The flow in this case is given by Poiseuille’s 
formula! 


Pi— +0.54 pum’, (1) 


where ~:— p2 is the pressure difference over the 
tube length L, ro the radius of the tube, p and 7 
density and viscosity of the gas, “», the maximum 
velocity in the center; at other places we have 


(2) 
with the average 
(2’) 


In a case such as the one treated by Rice,’ 
where f2=0.5 mm=660 dynes/cm’, L =30 
cm, 7o=0.5 cm, p=0.85X10-* (ethane, 0°C), 
one finds 


ti=u,,/2. 


Um ~ 1.2 104 cm ‘sec. (3) 


? Handb. d. Exper. Physik, Vol. IV-4, pp. 41, 46, Leipzig, 
1932. 

?F. O. Rice and Sister M. Denise Dooley, J. Am. 
Chem. Soc. 55, Oct. (1933). 


Even then, however, the flow is not turbulent, 
because the Reynolds* number is only 


R=itrop/n=50 (3’) 
(due to the small density) while the limit for 
turbulent flow is above 1000. 

In another case,‘ where the pressure drop is 
0.003 mm per cm, u»=2500 cm/sec., R=6.3 
(pressure 0.3 mm). 


B. HEAat CONDUCTION; GENERAL EQUATION 
AND SLOW FLow 


The equation of heat conduction is (A heat 
conductivity, C’ specific heat per cc). 


10T &T oT 


Or? r Or dx? Ox 


This equation has been integrated if wu is in- 
dependent of 7.5 

We can write \=aC’n/p where a is of the 
order of magnitude 1 (2.5 for a monatomic 
gas). 

The solution given in the reference above, 
however, is useful only if 


B= UoroC’ /2X= Uorop, 2an=R/2a (5) 


is not large compared with 2.4. This condition is 
not fulfilled by the conditions of the experiment 


3 Reference 1, p. 9. 

4W. R. Johnston, Diss., Johns Hopkins Univ. (1933). 

5 F. Paneth and K. F. Herzfeld, Zeits. f. Elektrochemie 
37, 577 (1931). 
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of Rice and Dooley, but is fulfilled in the case 
of Rice and Johnston.‘ Here 8=6.3/5~1.3 and 
the expression 


B— (6?+2.4*?)!= (6) 
The temperature difference falls off according to 


Q-2/1-15, 


This means that the gas has reached the wall 
temperature within 1 percent throughout while 
traversing a distance less than 3.5 cm. If R/2a 
is small compared with 2.4, then a distance of 
only /=2r9 would be required to reach tem- 
perature equilibrium up to 1 percent (this error 
would also affect the temperature coefficient to 
1 percent only). 


C. CONDUCTION; RAapip FLow 


Here, however, with R=50, the velocity is far 
too high to attain equilibrium over the whole 
cross section. Accordingly it is essential to take 
into account the fact that according to (2) the 
gas moves more slowly near the walls than in 
the body of the gas. We shall find in the following 
that there is a thin layer near the wall which 
closely attains temperature equilibrium, while 
the temperature of the core of the gas stream 
does not rise very high. In Rice’s experiment the 
test mirror is only a short distance (1 to 2 cm) 
beyond the end of the furnace. We shall attempt 
to give an approximate description of the 
mechanism according to the subsequent calcu- 
lations: Close to the walls there will be a thin 
sheet of gas in which the number of radicals 
corresponds to the temperature of the wall; 
the core will be relatively cold and will have only 
a negligible number of radicals. Because of the 
form of Arrhenius’ equation, with the rapid 
decrease of activated particles with falling tem- 
perature, only the highest temperature counts. 
At the mirror then, there will be a diffusion of 
molecules and radicals from the hot “‘layer’’ 
and of molecules from the ‘“‘cold’’ core. The 
latter act only as dilutants, and the amount of 
dilution is nearly independent of the temperature 
of the furnace. Accordingly, the temperature 
coefficient of the number of radicals is very 
little affected by the fact that they come from 
only a small fraction of the gas, though the 
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absolute number (which does not interest us 
here) is, of course, greatly diminished. 

The calculation is made in the following 
manner.® If heat conduction proceeds only in a 
thin sheet of gas adjacent to the walls, the 
curvature of this sheet can be neglected; that is, 
the problem can be treated as the heat con- 
duction by flow along a plate. This case has been 
calculated approximately by Pohlhausen.’ 

Let y be the distance from the wall (=7ro—r). 
Then his Fig. 2 (reference 7, p. 120) shows that 
(when 7,=temperature of wall, 7 )=original 
temperature of gas, and L’ = length of furnace) 


= —\), y= —(). 
nL’ roL’ 


(at about 1 mm distance, (7 7.) ~8 
percent according to his figure). 

If the number of radicals in 1 cc at temperature 
T; is given by 


(8) 
we find at the distance y 


N= = 1e-Cvu = (9) 


E 


D. THE DiFFUSION PROBLEM 


where 


C=0.3 


(10) 


We will consider next the problem of the 
diffusion of the radicals to the wall after leaving 
the furnace. Again, the equation of diffusion is 
given by (D=coefficient of diffusion) 


10N &N oN 
(— 
Ox? Ox 


(11) 


or? Or 


We can simplify this equation by considering 
the important sheet as a plane (neglecting 
(1/r)(@N/dr)); furthermore, we can neglect the 
effect of lengthwise diffusion (Ddé?N/dx*) as 
compared with that of convection (u@dN/dx), so 
that (11) takes the form 
D @#N/dy? =u dN/dx. (11’) 
*It has not been possible to integrate (4) suitably, 
taking account of (2). 


7E. Pohlhausen, Zeits. f. angew. Math. u. Mech. 1, 
115 (1921). 
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In integrating this equation, we will neglect the 
variation of u near the wall, that is put u= 4, 
which means a decrease of diffusion to the wall. 
Then (11’) has the familiar form of the equation 
of one dimensional heat conduction with No 
instead of 7, and x instead of ¢. Introducing the 
abbreviation 
u/4Dx =j, 

where x is the distance of the mirror from the 
furnace, the solution becomes 

(12) 


This, however, would be the distribution of No 
radicals starting from y=0, x=0. As the initial 
distribution at x=0 (end of the furnace) is 
given by (9), we find ultimately that 


N=N,(j/m)! f dy’, (13) 
0 


The number of radicals reaching a cm? of the 
wall per second through diffusion is 


or 


(14) 
0 
This can be transformed into 


(14’) 


The numerical value of this expression de- 
pends on 


To 2 R Dx 
roL 
To 2x pD 


Cc? E 
—=0.18( 
4j RT, 71 
E 
-0.18( 
RT, 7; 


But the theory of gases in its corrected form 
gives for the diffusion of an impurity (radical) 
present in small quantities*® 


Dp/n = 30/5 


(15) 


(15’) 


8’ D. Enskog, Phys. Zeits. 12, 56, 533 (1911); Ark. f. 
Math. Astr. och. Fys. 16, No. 16 (1922). 


K. F. HERZFELD 


where b is a numerical coefficient of order one 
defined by 


b= (1+ 


(d is the radius of collision; the primed quantities 
refer to the radical, the unprimed to the mole- 
cule). 
Therefore, introducing the abbreviation W, 
(W =C/2j') (15) becomes 


This has a value of approximately 2 in the case 
in which we are interested. With any value 
larger than 1, we have® for 
awe 


(16) 


The temperature coefficient of (17) is practically 
that of N, (see (8)), as the factor (1/W?)(1 
—3/2W?+15/4W?) changes approximately as 
T-?-*, The limits for our formula are: 


R<K1000 
2Rro/L'>1 
W=2. 


(15’’) 


R7, 


laminar flow, Eq. (3’) 
to give a thin hot sheet, (7) 
(15”’) 


Two more remarks must be made: We have 
treated only the heat transfer to translational 
and rotational degrees of freedom, which are 
measured by X, but not the vibrational, which 
are really responsible for the breaking of the 
bond. Now if, as is probable, this latter transfer 
is slow’ that does not alter our calculations, so 
long as its temperature coefficient is small (or 
entirely absent). For example, if it took a number 
N” of collisions with sufficient energy for the 


®E. Jahnke and P. Emde, Funktionentafeln, Leipzig, 
p. 31, 1909. 

10 If the heat transfer to vibrational degrees of freedom 
should be about as fast as the one to translational degrees, 
it would not affect the calculation for the rapid flow while 
it would increase the necessary length of the furnace in 
the slow flow somewhat. 
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excitation of the vibrations, this would diminish 
the absolute number but would not change the 
temperature coefficient. Secondly, in our calcu- 
lations we have neglected the fact that due to 
the change of temperature 7 and p change in the 
furnace. This will again not modify the tempera- 
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ture coefficient sufficiently to be of any sig- 
nificance. 

Of course, the velocities u are so high (close 
to sound velocity), that the use of the Poiseuille 
formula might be debatable, but nothing better 
is available. 
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Ionization Time of Thyratrons 


L. B. SNoppy, General Electric Research Laboratory 
(Received June 27, 1933) 


The ionization time for four types of commercial 
Thyratron tubes has been determined for the following 
cases: (1) impulsive voltage applied to the anode with 
the grid biased positively or negatively, (2) impulsive 
voltage applied to the grid with the anode potential 
constant, and (3) impulsive voltage applied to grid and 
anode simultaneously. The principal factors in fixing the 
ionization time are (1) the tube design; (2) the mercury 
vapor pressure (the time lag in starting is increased and 
made more erratic by decreasing the vapor pressure); 
(3) the grid bias (increasing the bias negatively increases 
the time lag); (4) the steepness of the wave front applied 


to the system. The lags become less as the rate of building 
up of voltage increases. This is due partly to the higher 
voltages reached and partly to the greater induced grid 
voltages when the impulse is applied to the anode. This 
effect is computed for two kinds of applied voltage waves. 
The time between grid and anode breakdowns was meas- 
ured for an FG-41 and an FG-17. The grid and anode 
currents during this interval were computed. For the 
particular FG-41 used, tanode = 3tgria. A table of ionization 
times under these various conditions is given for four 
Thyratron tubes (FG-17, 33, 41, 57). 


N all gaseous discharge devices initially in an 

unconducting state there is a definite time lag 
between the application of potential to the 
device and the attainment of maximum con- 
duction. During this time the gas path is ionized 
and the necessary emission built up at the 
cathode. In tubes containing emitting cathodes, 
such as the Thyratron tube, this delay allows 
the production and movement of positive ions to 
take place in sufficient numbers to neutralize the 
negative space charge at the cathode. A measure- 
ment of this lag in various types of Thyratron 
tubes was made some time ago and a brief report 
published.! In view of the increasing use of 
Thyratron tubes in scientific investigations it 
was thought worth while to give the detailed 
results. 

The time lag of Thyratron tubes has also been 
measured by Nottingham? and by Koch.’ The 
method used by Nottingham did not allow the 
actual delays to be measured, as intervals shorter 


1A. W. Hull and L. B. Snoddy, Phys. Rev. 37, 1691A 
(1931). 

2W. Nottingham, J. Frank. Inst. 211, 751 (1931). The 
long apparent time lags (of the order of 1000 10~* sec.) 
observed by Nottingham in his first paper on Thyratrons 
(J. Frank. Inst. 211, 271 (1931)) were due to an incorrect 
method of measurement. 

*W. Koch, Phys. Zeits. 33, 934 (1932). 


than about 20X10~* sec. could not be measured 
accurately. However he was able to determine 
that the actual lags were probably considerably 
less than this amount. The results of Koch are 
very accurate and of high time-resolving power. 
He obtained time lags ranging from about 
1.5X10-7 sec. up to a number of microseconds 
(10 to 20), the exact value depending upon the 
experimental conditions. However, they refer to 
specia! tubes and the method of applying the 
potential is not that ordinarily used with the 
Thyratron tube. In his work an impulsive po- 
tential was applied to the anode with the grid 
biased by different amounts from a constant 
voltage source. In the usual Thyratron tube 
circuits it is customary to apply the impulse to 
the grid with the anode kept at constant po- 
tential. The lag in the two cases may be quite 
different. Both of these methods were used in the 
present work. 

The impulsive potential was applied by a 
condenser circuit with a spark gap as the con- 
necting switch and the lags were determined by 
a cathode-ray oscillograph of the high voltage 
Dufour type (G.E. HC-2). 

The circuit shown in Fig. 1 was used for the 
determinations in which a positive potential was 
applied to the anode. The condenser C; was used 
as the source of potential. In all this work it was 
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Fic. 1. Circuit used to apply an impulsive potential to 
the anode of the Thyratron with the grid biased by the 
battery V, through the resistance R. 


i 


charged to 27 kv. The wave front was varied by 
changing Ri, Re, R; and C2. The rate of dis- 
charge through the tube after breakdown was 
limited by Ry. The grid was biased by the battery 
V, through the resistance R,. 

An oscillogram, typical of this arrangement, is 
given in Fig. 2. The kind of voltage wave 


6 & 
Kio” Secs. 


Fic. 2. Oscillogram showing the type of breakdown 
obtained with the circuit of Fig. 1. The applied voltage 
wave is shown at a and breakdown on this wave at b. 
Thyratron FG-17 with V,=—6.9 volts and R,=1000 
ohms. Condensed mercury temperature of 35°C. 


applied to the tube is shown at a, and the break- 
down on this wave at b. The tube used was an 
FG-17 with V, = —6.9 volts and R, = 1000 ohms. 

The circuit used to apply the impulsive po- 
tential to the grid with the anode at constant 
voltage is shown in Fig. 3. The grid was usually 
kept at —125 volts before the start of the im- 
pulse. With the oscillograph connected as shown, 
the time lag in grid breakdown could be de- 
termined but no indication of the anode break- 
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Fic. 3. Circuit for applying an impulse to the grid with 
the anode at constant potential. 


down would be obtained. Fig. 4a is an oscillogram 
typical of this method. The grid discharge starts 
at t, and is complete at fa. The tube is an FG-41 
with an anode potential of 125 volts. 

The total time of breakdown was measured by 
introducing a resistance of 1000 ohms at r. The 
current through this resistance raises the po- 
tential of the cathode and consequently that of 
the grid. In this way the two stages in the dis- 


> 


Fic. 4. a. Starting of grid discharge with the circuit at 
Fig. 3. The breakdown starts at ¢, and is complete at f2. 
Grid initially at —125 volts. Anode potential +125 volts. 
b. Grid voltage after addition of 1000 ohms at r. Grid 
discharge complete at t,—anode breakdown at t;. FG-41 
with mercury temperature at 40°C. 
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charge can be distinguished. This is illustrated in 
Fig. 4b. The grid discharge starts at ¢, and is 
complete at f2; the anode breakdown occurs at 
t;. An FG-41 was used for this picture with an 
anode potential of 870 volts and a mercury tem- 
perature of 40°C. The grid voltage rose to 220 
volts before breakdown. The total time of break- 
down is 4.9 10~* sec. 

This method was checked with the arrange- 
ment of Fig. 5. The spark gap G is coupled to 
the deflection plate circuit of the oscillograph 


OSC DEFLECTION 5000“ 
PLATES 
c 
2.5K10% 2 
- 
R 
+G- R Re 
218, > 
-OS 
= 
125 
mid vours 


Fic. 5. Circuit arranged to determine total time lag 
when impulse is applied to grid. Arc over of G registered 
on oscillograph plates by means of the small coupling 
capacity 


by the capacity C; (about 20uuf) and a 2.5 10° 
ohm resistance, as shown. Because of the in- 
ductance of the connecting wire P, sufficient 
voltage appears across the deflection plates to 
register the breakdown of G. With the wave 
front applied to the grid known, the total time 
lag can easily be determined. 

This is shown in Fig. 6 for an FG-41 with a 
constant applied anode potential of 125 volts 
and a mercury temperature of 35°C. The break- 
down of G is registered at ¢,; the discharge of the 
Thyratron tube starts at fg and is completed 
at f;. The interval t2—¢, is made up of two parts: 
(1) the time required to change the potential of 
the grid and its parallel 0.05uf capacity from 
— 125 volts to the breakdown grid potential; and 
(2) the time between grid and anode breakdowns. 
With the breakdown potential of the grid known 
for this particular wave the time between grid 
and anode discharge can be determined. 

The results of the tests are summarized in 
Table I. The voltages given are in every case the 
maximum values. In some of the work the grid 
and anode were connected together and the 


SNODDY 


Volts 


Rie secs. 


Fic. 6. Oscillogram of anode potential when using circuit 
of Fig. 5. Arc over of G at t;. Discharge starts at ts and is 
completed at ¢;. FG-41 with mercury temperature of 35°C. 


breakdown time determined for the combination. 
This is indicated in the table by equal grid and 
plate voltages. A grid bias and resistance (V, 
and R,) were used only in the cases in which an 
impulsive potential was applied to the anode. 
The negative values of V, represent the minimum 
grid potential at which consistent starting could 
be obtained with the given voltage impulses. 

By a study of the table the most important 
factor in determining the time lag for any par- 
ticular tube is seen to be the mercury vapor 
pressure. A decrease in pressure not only 
lengthens the starting time but also causes it to 
become more erratic. 

A negative bias also increases the lag as would 
be expected. On the other hand, the delay with a 
positive bias does not appear to be very sensitive 
to changes of several volts in the value of V,. 
With very steep waves the value of ’, seems to 
have very little to do with the lag, provided it is 
not negative enough to prevent starting. The 
reason for this is probably to be found in the 
finite inter-electrode capacities of the tubes. 


Fic. 7. Equivalent tube circuit. 


In Fig. 7 the equivalent tube circuit for an 
impulse applied to the anode is drawn with C; 
representing the anode-grid capacity, C2 the grid- 
cathode capacity and R the grid resistance. 
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TABLE I. Summary of tests on time lag of Thyratron tubes. 
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Time 
lag grid 
Time to plate 
Grid lag break- 
Tube Hg voltage - Anode grid down 
code pressure maxi- voltage (10-* (10-* Total time lag 
No. microns (volts) (ohms) mum maximum sec.) sec.) (10~* sec.) Remarks 
FG-17 6 +9.3 1000 462<+—+396 3.3<—+>2.8 Linear wave. AV/At 
= 1.65 X 10° volts/sec. 
6 +9.3 1000 924<—+550 1.2<—+0.7 AV/At=1.85 
X 10° volts/sec. 
15.5 +6.2 1000 176<—>110 1.3<—+1.0 1.65 X 108 volts/sec. 
4.1 +5.3 1000 77 4.0 Flat wave 
4.1 +2.4 1000 77 3.25 max. 
4.1 —6.9 1000 77 4.2<—+6.9 in 6.7 X1077 sec. 
8.8 44 264 i 0.95 0.95 Fig. 3. r= 10002 
8.8 66 66 0.83<—+0.95 Grid and anode 
connected 
FG-33 20.5 +6.25 1000 99 1.5 Flat wave 
20.5 +3.13 1000 99 1.8 
20.5 —1.6 1000 99 45 
20.5 88 88 1.1 Grid and plate 
connected 
4.1 +6.3 1000 99 6.4 Flat wave 
FG-57 20.5 —3.15 1000 99 3.1<——>2.7 Flat wave 
20.5 88 88 1.3 
4.1 99 99 12<—+>5.8 
FG-41 6 220 858(const.) 2.1 2.8<—-3 4.9<—-5.1 Fig. 3. r= 10002 
6 220 462 2.1 3.4 $5 Fig. 11 
6 176 125(const.) (3.2) 5 Fig. 5 
6 66 66 grid and plate 9.3<——+>8.6 8 X 10° volts/sec. 
6 154 154 connected 4.5 8 X10? volts/sec. 
4.1 176 125(const.) 2.2 5.4 8X 107 volts/sec. 
4.5 88 125(const.) a7 29 32.7 Flat wave 


With a voltage E = kt applied to this combination 
the potential of C2 is given by 


(1) 


This is plotted in Fig. 8 for an FG-17 and an 
FG-41. The constants used are given below: 


Tube Ci C2 R 
FG-17 4.7107" f 10.0 f 1000 ohms 
FG-41 1000 ohms 


The two rates of voltage increase given in 
Table I were used: k=1.65 10° volts/sec. and 
1.85 10° volts/sec. For an applied voltage of 


the form E = E,,4x(1—e~*) the potential of Ce is 
given by: 


EmaxaR Ci 


This is shown in Fig. 9 for the same Thyratron 
tubes and with two values of a: 0.25107 
sec.-' and 0.25 sec.-!. 


It is obvious that this induced grid voltage 
may be amply sufficient to start the Thyratron 
tube provided the grid is not initially negative 
enough to compensate for this effect. This is 
much more likely to be effective at high tem- 
peratures with the short time lags, especially so 
in the second case (Fig. 9) where the induced 
voltage is of relatively short duration. 

This difficulty may be overcome by connecting 
an additional capacity with short lengths of wire 
between the grid and cathode, as indicated by 
the dotted lines in Fig. 1. With an added capacity 
of 0.01yuf the rate of increase of Vc, for condition 
d of Fig. 8 is changed to that shown at e. The 
change in grid voltage is reduced to a negligible 
amount by this addition. For a potential wave of 
the form E= Eyax(1—e~**) the effect would be 
even more pronounced. 

It is evident that increasing the resistance in 
the grid circuit would raise the induced grid 
voltage with this type of impulsive starting 
(impulse applied to anode) and so decrease the 
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Fic. 8. Graph of induced grid voltage with E=kt 
applied to anode. e. 0.01 uf added capacity between grid 
and cathode. Otherwise same as d. 
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Fic. 9. Induced grid voltage with E=Emax(1—e“) 
applied to anode. Emax = 100 volts. 
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10. Induced grid voltage with E=kt and R of Fig. 7 
replaced by an inductance L=4X 1077 henries. 
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lag at high temperatures. This agrees with the 
work of Koch who obtained decreasing lags as 
he increased the grid resistance. 

When R of Fig. 7 is replaced by an inductance 
L, the grid voltage, for E=kt, is given by: 


sin (3) 


This is plotted in Fig. 10 for an FG-17 and 
for an FG-41 with k=1.85X10° volts/sec. and 
L=4X10-" henries. From this it is seen that 
the inductance of any ordinary loop of con- 
necting wire in the grid system would not have 
an appreciable effect upon the starting time 
wave fronts of this steepness. 

The current flowing in the grid and anode 
circuits between the times of grid and anode 
breakdowns (te—f; of Fig. 4b) can be computed 
from the circuit constants and the grid voltage 
measured on the oscillogram. The constants used 
for Fig. 4b are as follows: 


R, = 34,000 ohms C;=0.05uf 
R2= 14,000 ohms R= 10° ohms 
R;=5000 ohms r= 1000 ohms 
C2= 0.004 uf E=870 volts 


Since the time constant RC; is large, the influence 
of E may be neglected during the interval 
t3—te. The existence of a definite current to the 
anode during the time f;—/2 was established by 
changing C; from 0.05uf to 0.0008uf. An oscillo- 
gram with this value is given in Fig. 11. Although 
C3 was charged in both cases to the same voltage 
(870 volts), the maximum appearing at & in 
Fig. 11 is only 462 volts, while in Fig. 4b it is 
858 volts. The latter agrees within the experi- 
mental error with E. The decrease in voltage for 
the case of Fig. 11 is evidently caused by the 
anode current flowing during the interval ¢;—f2, 
discharging the condenser C; (R is too high to 
allow appreciable current from the battery £). 
The amount of charge lost by C; is 3.21077 
coulombs. This quantity corresponds to a change 
in potential of only 6 volts when C; is 0.05uf, 
which would not be detected by the oscillograph. 

The grid voltage from fz to fs; can be repre- 
sented by V,=a+Kt, where a and K can be 
measured from the oscillogram. Assuming that 
the charging current to C2 is constant during 
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Fic. 11. Oscillogram of grid voltage for circuit of Fig. 3 
with r= 1000 ohms. C;=0.0008 uf. 


this time, the net grid current is easily computed, 
with the result shown in Fig. 12a. 

The anode current may be calculated by two 
methods: (1), by assuming that the grid voltage 
measured on the oscillogram is given by V, 
(measured) =(20 volt grid-cathode drop)+i, 
1000+7ZanoaeX 1000; (2), by assuming that 
Zanode=Atgria Where A is a constant of pro- 
portionality. Since the total charge transferred is 
known, the plate current can be obtained from 
an integration of the grid current time curve. 
The quotient of 3.210-7 by the result of this 
integration gives A. For this particular tube 
Zanode=3 tgria- The results of the two methods 
are identical within the experimental error in- 
volved in reading the oscillogram. The anode 
current is plotted in Fig. 12b. The actual ratio 
between the two currents probably depends 
primarily upon the tube geometry. An FG-17 in 
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Fic. 12. a. Net current to grid for the interval t3;—t, 
computed from the oscillogram of Fig. 11. b. Anode 
current during same interval. 


the same circuit shows much less net grid current 
during the interval between grid and anode 
breakdowns, while the anode current is con- 
siderably larger. 

The lag between grid and anode breakdowns 
averages 0.53X10-* sec. less with C;=0.05uf 
than with C;=0.0008uf. This decrease is prob- 
ably due to the constant high plate voltage pre- 
vailing in the first case. 

The net current to the grid at the point 4, 
(Fig. 11) may also be computed. The circuit 
constants, however, were not measured with 
sufficient accuracy to allow anything more than 
an approximate estimate of this current to be 
made. All that can be said is that it does not 
exceed 8X10-* amperes, and it may be much 
less than this. 
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Ball Lightning 


J. C. Jensen, Department of Physics, Nebraska Wesleyan University, Lincoln, Nebraska 
(Received July 21, 1933) 


The size, appearance and physical properties of ball 
lightning has been subject to much controversy and but 
few accurate data are available. Description is given of a 
series of five photographs showing luminous masses which 


were of unusual size and which remained visible in changing 
form for about 3 minutes. Suggestions for further investi- 
gations are given. 


NE of the most difficult phenomena in 
nature to study is that known as ball 
lightning. It is so rare that comparatively few 
persons have ever seen it, and these observers 
have been so startled that they cannot give a 
clear description of what they saw. Even a 
trained scientist finds extreme difficulty in 
making any worth-while tests or measurements 
on a ball of blue fire as large as two fists which 
rolls along the floor for three or four seconds and 
then collapses with a noise resembling a big 
firecracker, leaving behind nothing more tangible 
than the odor of ozone. Yet these are the salient 
points in the descriptions reported recently by 
Talman,' Reynolds,? Mathias,* Marchant* and 
Craigie.° 
During the thunderstorm seasons of 1929, 
1930 and 1931, the writer was engaged in a 
study of the electrical field changes caused by 
lightning discharges with special reference to 
relation of the direction of branching of the 
flash to the sign of the field change.® This in- 
volved photographing a large number of dis- 
charges at night, most of them at close range as 
the storm center approached or came directly 
overhead. The path of the storm centers was 
traced by means of a battery of microbarographs. 


1C. F. Talman, Amer. Mercury 26, 69 (1932). 

2W. C. Reynolds, Nature 126, 413 (1930); ibid. 128, 
584 (1931). 

3E. Mathias, Comptes Rendus 194, 413 (1932); 194, 
2257 (1933); 196, 654 (1933). 

4E. W. Marchant, Nature 126, 128 (1930). 

5H. Craigie, Science 42, 344 (1930). 

6]. C. Jensen, Monthly Weather Rev. 58, 115 (1930); 
Phys. Rev. 40, 1013 (1932). 


During the early evening of August 30, 1930, 
after a sultry day characterized by local thunder- 
storms in the surrounding territory to the south, 
a severe storm of the “‘line-squall’’ type de- 
veloped in the west. The wind had been scarcely 
perceptible from the southeast but changed to 
the southwest and increased to twenty or thirty 
miles per hour between 9:35 and 9:45. Rain, 
accompanied by high northwest winds, began 
at 9:53. The barograph record shown in Fig. 1 
indicates an increase of 0.28 cm (0.11 inch) in 
pressure in the front of the storm, followed by a 


Fic. 1. Microbarograph record, storm of August 30, 1930 


second rise of 0.12 cm (0.05 inch) when the 
active part of the storm, with its rapid con- 
vection currents, passed over the instrument. 
The time on this record is about twenty minutes 
slow as the first rise began at 9:35. 

For making the photographs two cameras 
were placed in a fourth-story window above the 
trees and commanding a clear view of the western 
horizon. One of these is a 5X7 inch Graflex 
with a Bausch & Lomb Tessar lens rated at 
f 4.5, the other an Eastman 1A Kodak with 
f 6.3 lens. 
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BALL LIGHTNING 


Fic. 2. First 5X7 inch picture, enlarged. Fic. 3. First kodak picture, enlarged. Simultaneous with Fig. 2. 
Fic. 4. Second Graflex picture, enlarged. Fic. 5. Third 5 X7 inch picture, enlarged. 


The cold air rushing ahead of the nimbus 9:40 p.m. In the wake of one of these flashes 
cloud was filled with a swirling mass of dust but which, incidentally, carried negative electricity 
nevertheless brilliant flashes were seen descend- from cloud to ground, there appeared a shape- 
ing in rapid succession from the cloud to the less mass of lavender color which seemed to 
earth when the first films were exposed at float slowly downwards. The writer was so 
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occupied with the details of the photographic 
routine which required the operation of shutters, 
the use of a stopwatch for timing the lightning- 
thunder interval, entering data in the notebook, 
and telephoning to the assistant who had charge 
of the electrical recording apparatus in the 
laboratory, that there was little opportunity for 
close observation of this beautiful but unexpected 
phenomenon. The rose-colored mass seemed most 
brilliant near the ground and gave the im- 
pression of a gigantic pyrotechnic display. Two 
or three of the globular structures seemed to 
roll along a pair of 2300 volt power lines for 100 
feet or more, then bounced down on the ground 
and disappeared with a loud report. 

The first five of the fifty films exposed that 
evening cover an interval of about three minutes 
and all show the fire-balls in one of their stages. 
Fig. 2 is an enlargement of the first Graflex 
picture at 9:40 p.m. The one caught at the same 
instant with the Kodak, Fig. 3, resembles it 
closely, the lightning streamers in the rear being 
identical. Fig. 4 is an enlarged print from the 
second Graflex negative. It shows the principal 
descending ball in its most luminous and con- 
centrated form. Although the shutter was open 
for only a few seconds, two brilliant lightning 
discharges appear in the background, only one 
of which is included in the picture. That one, 
with streamers extending downward, gave a 
field change which indicated that negative elec- 
tricity had escaped from the cloud to earth. The 
Kodak shutter was closed while Fig. 4 was 
being taken but both cameras caught the final 
stage as shown in the enlargement of Fig. 5. 
The Kodak lens was left open longer and shows 
less detail, together with blurring, due to motion 
of the illuminated balls. 

As mentioned above, some of the fire-balls 
bounced off a power-line which was about 600 
meters (2000 feet) from the observer. From the 
well-known optical formula relating the size of 
object and image to object distance and focal 
length of the lens, the diameter of the first upper 
ball in Fig. 2 is shown to be 8.5 meters (28 ft) 
and that in Fig. 4, 12.8 meters (42 ft). Their 
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height above the horizon was 28 meters (92 ft). 
It will be objected that these dimensions are so 
much larger than the conventional size that they 
must relate to different phenomena. However, 
Toepler’ states the ‘‘masses of light’’ may vary 
in size from a hen’s egg to 10 meters in diameter. 
During the course of our researches balls of 
similar magnitude were seen on another occasion 
at a distance of two and one-half miles but no 
photographic record was obtained. In this case 
also, the balls collapsed with a sharp, loud re- 
port. 

A tornado which occurred on the evening of 
July 9, 1932, near Rock Rapids, lowa, gave 
evidence of a closely related type of luminous 
display according to the report of Mr. George 
Raveling, U. S. Weather Bureau observer. From 
the sides of the boiling, dust-laden cloud a 
fiery stream poured out like water through a 
sieve, breaking into spheres of irregular shape as 
they descended. No streak lightning of the usual 
type was observed and no noise attended the 
fire-balls other than the usual roar of the storm. 

‘While it is not the purpose of this article to 
propose a scientific explanation for the formation 
of the phenomenon described above, attention 
is called to the fact that ball lightning occurs 
most frequently in connection with dust from a 
chimney or fire-place, or in the dust-laden air 
of a wind-squall or a tornado. The interesting 
laboratory experiments of Cawood and Patter- 
son,* in which aerosols 20 cm in diameter and 
containing many of the characteristics of ball 
lightning have been produced are suggestive of 


’ further work along similar lines. An extension of 


the theoretical studies of Ollendorfi® on the 
pressure, temperature and conductivity relations 
in the luminous column of a lightning discharge 
would also be helpful in arriving at a better 
understanding of the physical factors involved. 


7Max Toepler, Mitteil. 
Isolatoren 25, 18 (1926). 

8 W. Cawood and H. S. Patterson, Nature 128, 150 
(1931). 

®F. Ollendorff, Archiv. f. Elektrotechnik 27, 169 (1933). 
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Velocity of Elastic Waves in Granite and Norite 


L. Don Leet, Harvard Seismograph Station* 
(Received June 17, 1933) 


The velocities of elastic waves were investigated in 
granite at Rockport, Massachusetts, and norite at Sud- 
bury, Ontario. The waves measured were generated by 
dynamite explosions, and recorded by portable seismo- 
graphs at distances up to 6260 feet in granite and 10,320 
feet in norite. The observed velocities in km/sec. with 
probable errors indicated by + were: 


Longitudinal Transverse Rayleigh 


Rockport granite 5.14+0.005 2.70+0.02 2.19+0.007 
Sudbury norite  6.22+0.003 3.49+0.006 2.79+0.01 


From these velocities, and the densities of specimens taken 
from the shooting locations, the following elastic constants 
may be computed. They are shown below, for comparison, 
with those computed by a similar procedure in an earlier 
report for Quincy, Massachusetts, granite. 


A quartz-rich zone in the norite was found to transmit the 
waves with a slower speed. The results are compared 
with statical measurements of Zisman on rocks from the 
same locations, a part of the same program of geophysical 
research at Harvard. Compressibility, measured directly, 
under hydrostatic pressures, on specimens exposed to the 
liquid transmitting the pressure is found to be the only 
statically determined constant which yields seismically 
effective values. The geological significance of granite 
velocities and elastic constants is discussed. It is concluded 
that a longitudinal-wave speed of 5.5 km/sec. cannot be 
regarded as unique for granite, or even necessarily sug- 
gestive of granite, as has been repeatedly proposed; and 
that there is as yet no seismological justification for the 
identification of any ‘‘layer’’ of the earth’s crust as either 
sedimentary or granitic. 


Quincy Rock- Aver- Sud- 
gran- port age bury 
ite granite granite norite 


Quincy Rock-  Aver- Sud- 
gran- port age bury 
ite granite granite norite 


Density 2.61 2.63 2.62 2.85 Young’s modulusX10~" 4.3 5.0 4.65 8.82 
(grams/cm*) (dynes/cm?) 

1.63 1.92 1.78 3.47 | CompressibilityX10" 22.8 22.8 228 15.6 
Bulk modulus X10 44 44 44 64 (cm*/dyne) 

(dynes/cm?) Poisson’s ratio 0.33 0.31 0.32 0.27 


I. LocaTION AND METHOD 


N a previous publication! some velocities in 

granite were reported. These included longi- 
tudinal waves in Rockport, Massachusetts, gran- 
ite. The present results represent an extension of 
that work at Rockport, together with a new 
investigation of the velocities in the norite of 
Sudbury, Ontario. Maps of the areas are shown 
in Fig. 3. The observing stations and shot 
locations shown there were mapped by plane 
table, at Rockport on a scale of 300 feet to the 
inch, and at Sudbury, 400 feet to the inch. 

The usual methods of seismic surveying were 
employed. The elastic waves were produced by 


*Paper No. 7, published under the auspices of the 
Committee on Geophysical Research and of the Division 
of Geological Sciences of Harvard University. 

1L. D. Leet and W. M. Ewing, Physics 2. 160 (1932). 


the explosion of charges of dynamite. Their 
arrival at a point at any desired distance from 
the explosion was recorded by a seismograph 
placed at that point. The distance travelled by 
the waves and the time required constituted the 
observed data. Distances were scaled from the 
originals of the maps reproduced in Fig. 3. 


II. INSTRUMENTS 


The seismographs used at Rockport have been 
described.' At Sudbury, electromechanical-trans- 
ducer seismographs, known commercially as 
geophones, were used. They, together with the 
necessary oscillographs, as well as radios for 
communication, were the generous gift of Dr. E. 
L. DeGolyer of the Geophysical Research Cor- 
poration. Vacuum-tube amplifiers were inter- 
posed between geophones and oscillographs. 
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Fic. 1. Specimens of Rockport records. 
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ELASTIC WAVES IN GRANITE AND NORITE 


The explosion of the dynamite was caused to 
interrupt a radio signal which was being recorded 
galvanometrically, thus registering the instant of 
the blast. Timing was effected by a tuning fork 
controlling the interruption of a light beam in 
the usual way. 


III. Rockport GRANITE RESULTS 


Specimen records from Rockport are repro- 
duced in Fig. 1 and the time-distance graph is 
shown in Fig. 2. Velocities were obtained by 
least square solutions. 


TIME-DI STANCE DATA 
FO 
ROCKPORT GRANITE 
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The recording distances of the previous work 
were extended to 6260 feet, leading to a revised 
figure for the longitudinal velocity of 16,840+20 
ft./sec. or 5.14+0.005 km/sec. 

The three-component seismograph, because of 
its low sensitivity, was not practicable for use 
throughout the investigation. It was, however, 
used at three different locations to give key 
records for the interpretation of the vertical- 
component records obtained elsewhere. These 
are shown as Nos. 3, 4 and 5 in Fig. 1. Of these, 
No. 4 shows the three principal phases practically 


a ° 
1 t 1 11 t 
23 2 6 10 " 12 RECORDS IN FIGURE 
1000 2000 3000 4000 3000 ahes 
DISTANCE iN FEET 
Fic. 2. 


undisturbed by natural frequencies of the re- 
cording location, which confused many of the 
records. The component writing the top trace 
was oriented to record only horizontal displace- 
ments in the plane perpendicular to the line 
joining shot and recorder, that is, transverse to 
the propagation direction. The middle trace rep- 
resents vertical displacements, and the bottom, 
longitudinal. 

It may be noted that the first arrival was a 
compression recorded on the longitudinal and 
vertical components. The beginning of motion on 
the transverse component was accompanied by 
some new motion on the vertical component. 
This indicated the feasibility of using only a 


vertical component seismograph even for the 
transverse wave. The transverse wave was 
followed, closely at the short distances used, by 
a typical Rayleigh wave in which the ground 
moved, with reference to the shot point, away- 
up-toward-down, the retrograde sense found to 
be characteristic of Rayleigh waves from earth- 
quakes.” 

By using these three-component records 
as keys to the interpretation of the others, 
the transverse and Rayleigh waves were plot- 


2L. D. Leet, Empirical Investigation of Surface-Waves 
Generated by Distant Earthquakes, Publ. of the Dominion 
Observatory 7, No. 6, 267-322, Ottawa (1931). 
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ELASTIC WAVES IN GRANITE AND NORITE 


ted. The velocities were, respectively, 8840+ 60 
ft./sec., or 2.70+0.02 km/sec., and 7190+20 
ft./sec., or 2.19+0.007 km/sec. 

If it is assumed that the transmitting medium 
obeys Hooke’s law, relations between wave 
velocities and elastic constants of the medium 
may be expressed by the following equations:* 


(1) 
Vr=(u/p)}, (2) 
k=r+2p/3, (3) 
B=1/k, (4) 
(5) 
(6) 


where V, is the longitudinal velocity, V7 the 
transverse velocity, k the bulk modulus, 8 the 
cubical compressibility, » the rigidity, E Young’s 
modulus, o Poisson’s ratio, and p the density. 

The values of these constants already reported 
for Quincy granite! together with those computed 
for Rockport granite are listed in Table I. 


TABLE I. Elastic constants of granite and norite. 


B 
E cm?/dyne 


k 

p (<10-" dynes/cm?) o 
Granite 2.61 1.63 33 44 4.3 22.8 0.333 
Quincy 
Granite 2.63 1.92 3.1 44 5.0 22.8 0.31 
Rockport 
Granite 2.62 1.78 3.2 44 4.7 22.8 0.32 
Average 
Norite 2.85 347 4.1 64 8.8 15.6 0.27 
Sudbury 


Several of the records of Fig. 1 show well- 
marked phases after the Rayleigh waves. In 
many cases, in fact, there appears to be an 
exact duplication of the first part of the record, 
as though there had been a double shot. Shooters 
have reported apparently double and sometimes 
triple shocks at the shot point. The possibility 
that these may be traceable to peculiarities of 


5A. E. H. Love, A Treatise on the Mathematical Theory 
of Elasticity, Cambridge University Press, Fourth Edition, 
643+ xviii pages (1927). 
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the quarry is lessened by the appearance of a 
similar phenomenon at Sudbury as shown on 
the records in Fig. 4. These later phases at both 
Rockport and Sudbury are as yet unexplained, 
and further investigation of them is proposed. 


IV. SupBurRY NoritTE RESULTS 


Copies of specimen records from Sudbury are 
reproduced in Fig. 4, and the time-distance 
graph is shown in Fig. 5. The velocities obtained 
from these data were: 


Longitudinal = 20,400 +10 ft./sec. 
Transverse =11,440+20 ft./sec. 
Rayleigh = 9,160+30 ft./sec. 


or 6.22+0.003 km/sec. 
or 3.49+0.006 km/sec. 
or 2.79+0.01 km/sec. 


A low-velocity zone was encountered in the 
first 1200 feet of the profile. Through the 
courtesy of Dr. Arthur Yates, chief geologist of 
the International Nickel Company, at Copper 
Cliff, data on this situation were obtained from 
a suite of thin-sections taken along a 3700-foot 
line northwest from the Murray Mine. Examina- 
tion of these showed a sudden drop in the 
percentage of quartz at about the neighborhood 
of 1200 feet from the edge of the norite. From 
twenty percent, it dropped to five and less. It 
has been found by L. H. Adams in laboratory 
studies that ‘‘a decreasing silica content in rocks 
is accompanied by a steady decrease in com- 
pressibility,”* which confirms this finding of 
lower velocity with higher percentages of quartz. 

Elastic constants of the norite, computed from 
velocity data, are listed in Table I. Dr. D. S. 
Muzzey of the Harvard Research Laboratory of 
Physics has obtained a comparable value for E 
by using a magnetostriction method to set the 
specimen vibrating. The mean of his values on 
four Sudbury norite specimens is 9.110" 
dynes/cm?. 

As in the case of Rockport granite, there is a 
second large phase which is unexplained as yet. 
As nearly as its apparent beginning can be 
picked, it plots with a velocity of about 7550 
ft./sec. or 2.30 km/sec. 


4L. H. Adams, Gerlands Beitrage zur Geophysik 31, 
315-321 (1931); also in Papers from the Geophysical 
Laboratory, Carnegie Institution of Washington, No. 770. 
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Fic. 4. Specimens of Sudbury records. 
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V. ANALYSES AND STRUCTURAL RELATIONS 


Chemical and mineralogical analyses of Rock- 
port granite have been published,’ as have 
chemical analyses of Sudbury norite.® 

The following is a mineralogical analysis of 
the norite, made by Dr. H. J. Fraser when a 
member of the Harvard University Department 
of Geology. Specimens were obtained through 
the courtesy of the International Nickel Com- 
pany. Volume percentages of the constituent 
minerals are given in Table II. 


TABLE II. Volume percentages of minerals—Sudbury norite. 


Core Core Core Core 
1 3 4 
Labradorite, AbysAns2 67 58 70 
Andesine, AbssAna,4 74 — 
Potash feldspar 1 1 2 4 
Hypersthene 
(25 percent MgSiO;) 15 8 18 10 
Augite 3 5 4 a 
Hornblende 2 5 8 3 
Biotite 1.5 6 1.5 3 
Quartz 3 8 8 6 
Magnetite, apatite, sericite, 
kaolin, carbonate 0.5 oa 0.5 1.5 


Barlow,® Coleman’ and Phemister* have pub- 
lished reports on the Sudbury area. Coleman’s 
report is accompanied by a useful large-scale 
map, and Phemister’s by a small-scale sketch of 
the general area. Both show the Murray and 
Elsie mines indicated on Fig. 3. 

Because of its association with ore-bearing 
zones, the norite has been carefully studied in 
this area, and its structural relations are well 
known in many places. In locating contacts of 
the norite, we were materially aided by advance 
maps embodying the results of a remapping of 


5 Coleman, Moore, Walker, 
Studies, No. 28, p. 18 (1929). 

6° A. E. Barlow, Origin, Geological Relations, and Compo- 
sition of the Nickel and Copper Deposits of the Sudbury 
Mining District, Geological Survey of Canada, Publ. No. 
961, Ottawa (1907). 

7A. P. Coleman, The Nickel Industry: With Special 
Reference to the Sudbury Region, Ontario, Canada Depart- 
ment of Mines, Mines Branch, Publ. No. 170, Ottawa 
(1913). 

’T. C. Phemister, Zgneous Rocks of Sudbury and their 
Relation to the Ore Deposits, 34th Annual Report of the 
Ontario Department of Mines, Vol. 34, Part VIII, Toronto 
(1925). 
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the area which has been carried on since 1930 
under the direction of the Geological Survey of 
Canada. For these, some excellent aerial photo- 
graphs of the zone of the velocity determina- 
tions, and valuable suggestions, I am greatly 
indebted to Dr. W. H. Collins, Director of the 
Survey. 

Through the further generous cooperation of 
Dr. Yates of the International Nickel Company, 
data from a number of bore-holes_ located 
approximately along our velocity profile were 
obtained. Fig. 6 is constructed from these, 
showing the dip of the bottom contact of the 
norite. 


VI. Discussion 


Attention was called in a previous report! to 
the great disparity between the elastic constants 
of Quincy granite computed from wave speeds 
observed in the field and those derived from 
laboratory experiments of Adams and Coker® at 
relatively low pressures. Since that time, under 
the direction of the Harvard Committee on 
Geophysical Research, Zisman has performed an 
elaborate series of experiments which clarify the 
situation. Published values of such elastic con- 
stants as Young’s modulus, Poisson’s ratio, and 
the compressibility, determined by statical meth- 
ods for similar rock types, have ranged between 
wide limits. For granite, Young’s modulus was 
reported by the Watertown Arsenal!® as ranging 
from 3X10" to 6.510" dynes/cm?; by Adams 
and Coker® from 4 to 6.3; by Nagaoka" from 
1.1 to 4.2; and by Bach” for German granites as 
2.4. The Watertown Arsenal reported Poisson’s 
ratio for granite between 0.07 and 0.25, while 
Adams and Coker observed values between 0.20 
and 0.26. Compressibility of granite at low 
pressures was reported by Adams and Coker 
ranging from 30.3 10-" cm?/dyne to 36.8; and 


°F, D. Adams and E. G. Coker, Carnegie Inst. Wash- 
ington, Publ. No. 46, 69 pages (1906). 

10 Watertown (Massachusetts) Arsenal Annual Report 
for 1895. 

1H. Nagaoka, Elastic Constants of Rocks and the 
Velocity of Seismic Waves, Publs. of the Earthquake 
Investigation Committee, No. 4, 47-67, Tokyo (1900). 

12C, Bach, Elastizitét und Festigkeit, Fourth Edition, 
p. 64 (1904). 


ELASTIC WAVES IN GRANITE AND NORITE 383 


SURFACE 


EL. East 


| 
28 
N R ' T 
4 
4 
, 
4 
4 
4 
_- [963 
os 
1080 1099 
“> 
4250 
aaa 
a 
SUDBURY, ONTARIO 
at SECTION N 26° W OF ELSIE EAST PIT 
_-} 1640 
- 
73 00 FEET 
A BASIC FLOW 
ROCK 
491 
1631 1840 
“9 
2018 
Fic, 6, 


by L. H. Adams and Williamson" at pressures 
around 2000 atmospheres as 19.8 to 20.1. 
Zisman,'* working at low pressures, found for 
all rock types studied that Young’s modulus and 
Poisson’s ratio determined in air-dried specimens 
varied greatly with the mean stress used. The 
same was true for the compressibility of speci- 
mens covered to exclude the liquid which trans- 
mitted the pressure.'° He pointed out, as a result, 
that these variations with pressure have doubt- 
less caused much of the apparent disagreement 
among the results of previous experimenters, 
though he extended the work to show that this 
has not been the entire source of trouble. 
Adams and Coker® worked for the most part 
with mean stresses of about 315 kg/cm*®. They 
measured the longitudinal compression and 


3L. H. Adams and E. D. Williamson, J. Frank. Inst. 
195, 475-529 (1923). Reprinted as No. 484 of Papers from 
the Geophysical Laboratory, Carnegie Inst. of Washington 
(1923). 

4 W. A. Zisman, Young's Modulus and Poisson's Ratio 
with Reference to Geophysical Applications, Proc. Nat. 
Acad. Sci. 19, 653 (1933). 

1% W. A. Zisman, Compressibility and Anisotropy of Rocks 
at and Near the Earth's Surface, Proc. Nat. Acad. Sci. 19, 
666 (1933). 


lateral extension under a unidirectional thrust, 
from which the compressibility was computed. 
Zisman, in addition to determining Young's 
modulus and Poisson’s ratio by a similar method, 
measured compressibility directly under hydro- 
static pressures. The only point at which approxi- 
mately direct comparison is possible shows the 
results of these similar statical methods on 
Young’s modulus to be somewhat variable. 
Zisman’s figure for Young’s modulus in Rockport 
granite, extrapolated from 202 kg/cm? mean 
stress to 315 kg/cm* is 7.2010" dynes/cm? 
whereas that of Adams and Coker for Quincy 
granite at the same stress was 4.64 for one 
specimen and 5.68 for another. Zisman, further, 
found that Young’s modulus for Quincy granite 
varied greatly according to the location from 
which it was sampled, as had been suggested.! 
A surface specimen, for example, showed 2.13 
X10" dynes/cm? at a mean stress of 11.2 kg/cm? 
whereas one from a depth of 100 feet in a working 
quarry showed 5.18 at the same stress, and one 
at 235 feet in the same quarry, 4.85. These 
differences decreased with increased stresses. 
Zisman established an important point re- 
garding indirect determinations of compressi- 
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bility."° For Rockport granite, he found the com- 
pressibility calculated from Young’s modulus and 
Poisson's ratio at 11.2 kg/cm? to be 66.2 10-" 


cm?/dyne, while the value observed directly’ 


under hydrostatic pressures was 84.0. At 56.0 
kg/cm? the values were 54.5 and 65.0, respec- 
tively. For Quincy granite, the figures were 70.6 
as against 86.0 at 11.2 kg/cm? and 47.5 as 
against 67.0 at 56.0 kg/cm?. 

In the direct determinations of compressibility 
under hydrostatic pressures, all rocks when 
covered were found to exhibit rapid changes of 
compressibility with pressure until fairly ‘well 
compacted (at pressures near 1000 kg/cm?). The 
values for Quincy and Rockport granite un- 
covered, on the other hand, agree well with 
figures of L. H. Adams and Williamson for 
pressures from 2000 to 10,000 atmospheres,” 
and with values computed from elastic-wave 
velocities. Table III compares these values. 
These figures are in excellent agreement as con- 
trasted with such numbers as 91.7 for covered 
Rockport granite at zero pressure. 


TABLE III. Mean compressibility of granites. 


Laboratory determinations Computed from 
by direct measurement on _ elastic-wave 


Stress uncovered specimens velocities 
kg/cm? cm?/dyne X 1078 
Adams and 
Zisman Williamson 
0 21.1 
22.8 
720 18.1 
2000 21.1 
10000 18.3 


Zisman has stated that the presence of pores 
and cracks explains anomalous behavior of rocks 
under pressure until such discontinuities are 
closed. It is not surprising, then, that a liquid in 
the openings causes the rock to possess constants 
found in the field, where the elastic-waves 
studied travel almost entirely through water- 
soaked materials. Accordingly, measurements of 
compressibility under hydrostatic pressures on 
liquid-soaked rocks appear to be the only avail- 
able statical values at low pressures which yield 


*W. A. Zisman, Comparison of the Statically and 
Seismologically Determined Elastic Constants of Rocks, Proc. 
Nat. Acad. Sci. 19, 680 (1933). 
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seismically effective data. There is still no con- 
clusive evidence for a fundamental difference 
between statically and dynamically determined 
elastic constants for a given rock under com- 
parable conditions. Also, it appears that the 
extremely valuable law of L. H. Adams and 
Williamson: ‘‘the compressibility of a fresh holo- 
crystalline rock is an additive function of the 
compressibility of its minerals,’’ needs not be 
qualified as valid only ‘“‘at moderately high 
pressures,’’ but may be regarded as true for all 
pressures if voids are either absent or liquid- 
soaked. 


Geological significance of granite velocities and 
elastic constants 


Accumulating data on elastic constants and 
wave speeds indicate that great care needs to be 
exercised in attempting to use them for the 
identification of buried rocks, particularly near 
the surface. 

Jeffreys adopted the symbol P, for longi- 
tudinal waves from near-earthquakes, which 
show speeds averaging approximately 5.5 km/ 
sec. In so doing, he stated, ‘‘there no longer seems 
to be any reason for diffidence in supposing these 
waves characteristic of granite.’’ Many investi- 
gators have followed him in this and perpetuated 
the term “granitic layer’’ which, geologically, is 
a decidedly unfortunate generalization. 

The basic criteria for identification of surface 
“‘layers”’ of the earth have been wave speeds and 
comparisons of computed elastic constants with 
those measured in the laboratory. Field investi- 
gations have shown now, however, that longi- 
tudinal elastic-waves have speeds of the order of 
5.0 to 5.2 km/sec. in granite near the surface. 
On the other hand, 5.6 km/sec. has been re- 
ported by Brockamp and Wéolcken"’ as the 
velocity of longitudinal waves in a basalt; 
Richards'* found velocities of 5.18, 5.36, 5.64, 
5.97, and others in that range for limestone. 
Similar velocities have been recognized for lime- 
stone in commercial seismic prospecting in the 
Gulf Coast of the United States. Preliminary 
results (obtained in connection with the present 


17 Brockamp and Wolcken, Zeits. f. Geophysik 5, 163- 
171 (1929). 
18 T. C. Richards, Proc. Phys. Soc. 45, 70-81 (1933). 
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investigation) of an incomplete profile in a slate 
in Pennsylvania indicate a longitudinal velocity 
of the order of 6 km/sec. Zisman’s results show 
a compressibility for Pennsylvania limestone 
exactly the same as that for one surface specimen 
of Quincy granite, and the compressibility of a 
quartzitic sandstone only about ten percent 
greater. Richards got 22.810~-" cm?/dyne for 
the compressibility of one of his limestones. 

The new facts are serving rather to complicate 
than to clarify the situation at present. The 
ambiguity of current lithological interpretations 
of velocities and seismically effective elastic 
constants is looming larger at every turn. It 
cannot be said that a longitudinal wave speed of 
5.5 km/sec. is unique for granite, or even 


necessarily suggestive of granite. There is as yet 
no seismological justification for the identitica- 
tion of any “‘layer”’ of the earth's crust as either 
sedimentary or granitic. 
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